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The wetting behavior of water on the solid surfaces is fundamental to various physical, chemical and 
biological processes. Conventionally, the surface with charges or charge dipoles is hydrophilic, whereas the 
non-polar surface is hydrophobic though some exceptions were recently reported. Using molecular 
dynamics simulations, we show that there is a critical length of the charge dipoles on the solid surface. The 
solid surface still exhibited hydrophobic behavior when the dipole length was less than the critical value, 
indicating that the water molecules on the solid surface seemed not "feel" attractive interactions from the 
charge dipoles on the solid surface. Those unexpected observations result from the collective interactions 
between the water molecules and charge dipoles on the solid surface, where the steric exclusion effect 
between water molecules greatly reduces the water-dipole interactions. Remarkably, the steric exclusion 
effect is also important for surfaces with charge dipole lengths greater than this critical length. 

A microscopic understanding of the wettability of solid surfaces is fundamental to various physical 1-20 , 
chemical 21 " 27 and biological processes 28 " 33 , such as microfluidics 34 ' 35 and sensing 36 ' 37 . Particularly, the 
hydrophobic effect plays a key role in protein folding 38 " 40 , nano-scale dewetting transition 41 " 43 and self- 
assembly of amphiphiles 44,45 . However, the physics of the hydrophilic or hydrophobic effects is not fully under- 
stood because the hydrophilicity or hydrophobicity of a surface is usually subject to various factors, including 
temperature 46 and morphology 47 . Conventionally, the surface with charge dipoles is hydrophilic, whereas the 
non-polar surface is hydrophobic. However, some exceptions were recently reported. Zangi and Berne 48,49 showed 
that a high value (e.g., 1.5 e) of electrolyte particles in solution can induce hydrophobic particles to be more 
hydrophobic. Giovambattista et al. 50 performed a series of simulations on polar hydroxylated silica (Si-OH) 
surfaces, finding that the surface is macroscopically hydrophobic, even though the dipole moment magnitude 
equaled to 41% that of the water molecule. The topography of the surface was also found to play a significant role 
in the wetting behavior of the protein surface or hydroxylated silica surface 51 ' 52 . Boron-nitrogen nanotubes with a 
partial charge of ~0.4e 53 have been shown to still exhibit hydrophobicity 54 . Moreover, in 2006, Li et al. 55 showed 
that van der Waals interactions are the dominant interactions between the superhydrophobic surface molecule 
CF 3 (CF2) x (CH2) y and the water molecules on the surface with a charge of — 0.2 e in CF 3 groups. Very recently, we 
showed that hydrophobic-like behavior occurs on extremely polar hexagonal surfaces due to the ordered water 
monolayer for very large partial charges 4 ' 56 , which was recently confirmed experimentally 57 . However, a relatively 
complete picture about the effect of the surface charges and charge dipoles on the wetting property is still lacking. 

Here, using molecular dynamics simulations, we show that the length of the charge dipoles on the solid surface 
plays a key role on the wetting behavior. Moreover, there is a critical length of the charge dipoles on the solid 
surface, below which the charge dipoles on the solid surface plays unexpectedly negligible role in the wetting 
property. The solid surface still exhibits hydrophobic behavior when the dipole length is less than the critical 
value, indicating that the water molecules on the solid surface seem not to "feel" attractive interactions from the 
charge dipoles on the solid surface, no matter how large the moment magnitude of the charge dipoles is. Those 
unexpected observations result from the collective interactions between the water molecules and charge dipoles 
on the solid surface, where the steric exclusion effect between water molecules prevents those hydrogen atoms of 
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water molecules from staying very close to the negative charge and 
those oxygen atoms of water molecules from staying very close to the 
positive charge, reducing the interactions between the water mole- 
cules and the charge dipoles. Interestingly, the steric exclusion effect 
was also important for surfaces with charge dipole lengths greater 
than this critical length. When the charge dipole lengths are greater 
than this critical length, it is found that the surfaces become more 
hydrophilic with the dipole lengths increase. We note that the bond 
lengths of most of common materials are less than the critical 
lengths, suggesting that the observation should be of general import- 
ance to the surface wetting properties. 

The geometry of the solid surface system is shown in Fig. l(a,b) 
together with snapshots of some water molecules. Positive and nega- 
tive charges of the same magnitude, q, were assigned to every bond- 
ing atom of the hexagonal arrangement, represented by green and 
blue triangles, respectively, which was similar to the graphene and 
boron- nitrogen monolayer. The surface was neutral. The length of 
charge dipoles was denoted as /. Initially, 924 water molecules in a 
rectangular shape were located on solid surface composed of 2376 
atoms in a box. Each system was simulated by molecular dynamics 
(MD) for 10 ns, and the data of the last 2 ns were collected for 
analysis. For all simulations, the water slab was in phase-level coex- 
istence with the vapor at 300 K. 

Results 

Generally, the existence of polarity on the surface enhances the 
hydrophilicity. Surprisingly, our simulation shows that the contact 
angles of the liquid water droplets on the solid surfaces are always 
close to, even larger than 90° when / < 0.202 nm, for all q value in the 
interval of [0, 1.0 e]. In Fig. 1(c), we show the contact angle discrep- 
ancy (9 0 —9 q ) values associated with the parameter / for different q, 



where 0 q and 0 0 are the contact angles of the liquid droplets on the 
surface with the charge of q and without any charge, respectively. For 
/ < 0.202 nm, 0 q is very close to 0 0 for all charges less than 1.0 e, 
indicating that the effect of the charge dipoles on the wetting prop- 
erty of the surface is negligible. From Fig. 1(c), one can observe that 
0 0 — 6 q increases as / increases from 0.202 nm. The larger the value of 
q, the faster 0 0 — 6 q increases. The increase of 0 0 — 0 q is very quick for 
q > 0.6 e, and it reaches its maximal value, which corresponds to the 
water completely wetting the solid surface (6 q = 0), for q = 1.0 e at 
/ > 0.242 nm, for q = 0.8 e at / > 0.262 nm, and for q = 0.6 e at / > 
0.402 nm. A longer dipole length is clearly required to achieve com- 
plete wetting of the solid surface for a smaller charge. However, we 
find that the water could not complete wet the solid surface for q < 
0.4 e, even if / is very large (0.482 nm). Interestingly, for q = 0.4 e, the 
contact angle increases significantly only from / ~ 0.3 nm. Thus, the 
wetting properties of the solid surface significantly depend on the 
length of the dipole. We note that, when / < 0.202 nm, the surface 
always shows hydrophobic (6 q ~ 90°) when the dipole moment 
magnitude is very large, even larger than 3 times of that of the water 
molecule. In order to better characterize the behavior, we define a 
critical length as the dipole length below which the contact angle 
discrepancy (0 0 —0 q ) is less than a critical value, i.e., 15°, for all q 
value in the interval of [0, 1 .0 e] , or in the interval of [0, 0.4 e] , denoted 
by l c and / c 04 , respectively. l c = 0.202 nm and l c 0A = 0.302 nm, 
which are shown by grey and light blue shadows in Fig. 1(c), respect- 
ively. We note the bond lengths of many common materials are 
smaller than l c = 0.202 nm and most of common materials are 
smaller than / < l c 0A = 0.302 nm (see Fig. Id and the Table SI in 
the Supporting Information). This also explains that boron-nitrogen 
nanotubes still exhibit hydrophobicity although there are partial 
charges of ~0.4e on the boron nitrogen atoms 53 . 
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Figure 1 | Introduction to the model system and its wetting behaviour. (a,b) Structure of the model surface with /= 0.142 nm and 0.262 nmatg= 1.0 e 
together with snapshots of some water molecule. The green and blue triangles represent the atoms on the solid surface with positive charges and negative 
charges, respectively. The oxygen and hydrogen atoms of the water molecules are shown by (bright or soft) red and white spheres. The oxygen atom of a 
water molecule with one OH bond pointing downward is shown by soft red. (c) Contact angle discrepancy (9 0 ~ 0 q ) with respect to the dipole length (/) for 
the charge (q) 0.2 e (black), 0.4 e (red), 0.6 e (blue), 0.8 e (green), and 1.0 e (magenta). 0 q and 0 o are the contact angles of the liquid droplets on the surface 
with a charge of q and without any charge, respectively. The saturation values of 0 0 — 0 q increases from —102° for / = 0.242 nm to —130° for / = 0.482 nm 
due to the modification of the lattice length, which correspond to complete wetting of the surface (6 q = 0). (d) Critical lengths for some typical surfaces of 
binary compound crystals. Here, the critical length is defined as the dipole length below which the contact angle discrepancy (9 0 —6 q ) is less than 15° for all 
g value in the interval of [0, 1.0 e] , or in the interval of [0, 0.4 e] , denoted by l c and l c ° A and shown by black and red columns, respectively. The first surface 
type is the structure shown in (a,b). The other surface types, increasing from 2 to 6, corresponds to the square, rhombic and rectangle lattice, (110) surface 
of a face-centered cube, and (110) surface of a body centered tetragonal, respectively (see Fig. SI in Supplementary Information). 
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Discussions 

In order to understand the mechanism underlying this behavior, we 
have calculated radial distribution functions (RDF) between water 
oxygen (O) and hydrogen (H) atoms and the surface positive and 
negative charged atoms, denoted as go+(r), gH+(r)> go-(r) and 
g H -(r) 9 respectively. Some typical results for q = 0.4 e are shown 
in Fig. 2. Interestingly, the RDF profiles are almost independent of 
the charge positivity or negativity; go+{r) ~g 0 -(r), and g H +(r) ~ 
gH- (r) for all rvalues when / is very small (0. 142 nm). The peaks of all 
the RDF profiles are at r = —0.4 nm. We have shown that the RDFs 
of the O and H for q = 0 are consistent with the RDFs for q = 1.0 e 
when / = 0.142 nm (see Fig. SI in the Supporting Information). In 
fact, we find that the RDF profiles are almost independent of the 
charge for q < 1 .0 e, suggesting that the dynamics of the surface water 
is independent of the quantity of the charge dipoles on the surface for 
/ = 0.142 nm. In Fig. 1(a) we show the snapshot of the water mole- 
cules in the layer with a thickness of 0.5 nm from the top view for / = 
0.142 nm at q = 1.0 e. These water molecules do not show any 
ordering structure, indicating that water molecules seem not "see" 
the surface charges even q = 1.0 e. The RDF profiles for / = 0.182 nm 
were quite similar to the RDF profiles for / = 0.142 nm. When / = 
0.222 nm, a clear but small peak appears at r = 0.20 nm for g H -(r). 
At / = 0.362 nm, clearer peaks appears at r ~ 0.18 nm for g H -(r) 
profiles and at r ~ 0.27 nm for go-(r) profiles. The two peaks 
become clearer for / = 0.402 nm. The distance between those two 
peaks is found to be about 0.09 nm, which is very close to the bond 
length of O-H in the SPC/E water model of 0.1 nm used here. Careful 
examination finds that most of the O atoms near the peak at r ~ 
0.27 nm chemically bond H atoms near the peak r ~ 0.18 nm. This 
explains the observation that the peak in the RDF profile of go-{r) 
moves towards left as / increases since the first peak of gn-(r) moves 
towards left. From the snapshot for / = 0.262 nm atq = 1.0 e shown 
in Fig. 1(b), one can observe that the water molecules close to the 
solid surface show clear structure ordering, i.e, the water molecules 
above the negatively charged atoms have one OH bond pointing 
downward to the surface, and the water molecules above the posi- 
tively charged atoms have neither OH bond pointing downward, 
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Figure 2 | Radial distribution functions (RDF) of water oxygen O and 
hydrogen H atoms with respect to the surface positive and negative 
charged atoms for the surfaces considered, respectively, for various dipole 
lengths / at q = 0.4 e. The black lines and green lines correspond to the 
RDF profiles for the oxygen (O) and hydrogen (H) to the positive charges, 
denoted by g Q+ (r) and g H+ (r), respectively. Blue lines and red lines 
correspond to the RDF profiles for the oxygen (O) and hydrogen (H) to the 
negative charges, denoted by g Q - (r) and g H - (r), respectively. 
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indicating that the oxygen atoms are nearest to the positively charged 
atoms on the surface. Thus, the water molecules on the surface with 
/ = 0.262 nm behaves qualitatively different from the water mole- 
cules on the surface with / = 0.142 nm, consistent with the RDF 
profiles. There is also a peak for the RDF profile of gn+(r) at r = 
0.32 nm. As / increases more, the locations of the peaks in the profile 
of go+{r) moves towards smaller values of r, whereas the profiles of 
gH+( r ) show different characteristics. All of these observations show 
that the charge dipoles on the sold surface have a more profound 
effect on the distribution of the water orientation on the solid surface 
when / is larger. Due to the small van de Waals radius of an H atom, H 
is much closer to the negative charges compared to the distance from 
the positive charges to O atoms. 

The interactions between water molecules and the solid surface are 
usually regarded as the key factor in determining the wetting prop- 
erty of the surface. We have computed the average electrostatic 
interaction energy between each single charge (average value from 
both cases of positive and negative charges) on solid surface and all 
the water molecules near the charge, with respect to / for various q, 
denoted as E dectro (see Fig. 3). E dectro is very small for / < l c (the the 
grey shadows in Fig. 3). The water on the solid surface seems not to 
"feel" the attractive interactions from the charge on the solid surface, 
contrary to intuition, but consistent with our simulation results that 
solid surfaces with charge dipoles still exhibit hydrophobic behavior 
for / < l c . As / increases from l c , the electrostatic interaction increases 
monotonically. We have also computed the electrostatic interaction 
energies from the RDFs of water oxygen O and hydrogen H atoms, 
denoted by E theory . Herein, the partial charges on an oxygen and 
hydrogen atoms are —0.8476 e and 0.4238 e, respectively, from the 
SPC/E water model. As shown in Fig. 3, E theory agrees with E dectro very 
well. 

Now we may capture the physics behind this unexpected behavior. 
The key is the collective effect of the interaction of the water mole- 
cules with the charge dipoles on the solid surface, where the steric 
exclusion effect between the water molecules may become of import- 
ance. When / is large enough, the negative charges on the solid 
surface attract the H atoms and the positive charges attract the O 
atoms so that some of the H atoms are very close to the negative 
charges and some of the O atoms are very close to the positive 
charges. When / decreases to a smaller value, the steric exclusion 
effect prevents those hydrogen atoms from staying very close to 
the negative charge and those oxygen atoms from staying very close 
to the positive charge. The electrostatic interactions between water 
molecules and surface charge dipoles decrease correspondingly. 
When / is small enough, it becomes almost impossible to arrange 
the positions and orientation of the water molecules so that the O 
atoms are close to the positive charges while the H atoms are close to 
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Figure 3 | Electrostatic interaction energy between water molecules and 
each charge on the solid surface with respect to the dipole length (I) for the 
charge (q) 1.0 e (magenta), 0.8 e (green), 0.6 e (blue), 0.4 e (red), and 0.2 e 
(black). The open circles are determined directly from numerical 
simulations (denoted by E dectro ) and the filled stars represent the 
electrostatic interactions energy calculated from the radial distribution 
function (denoted by E f /j eory ). 
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the negative charges. It is clear that this value of / should be compar- 
able or smaller than the size of the water so that the positions and 
orientations of the water molecules near the solid surface cannot be 
arranged according to the charges on the surfaces. Our numerical 
results have shown that the critical surface dipole length l c = 
0.202 nm and l c 0A — 0.302 nm are close to the diameter of a water 
molecule (—0.28 nm). Consequently, when / < / c , the atom densities 
are almost independent of the signs and the qualities of the charges 
and the water molecules on the solid surface seem not to "see" the 
charges on the solid surface. 

The existence of the critical length is universal. We have studied 
the wetting behavior of the surfaces of other five typical types. Three 
of them have 1 : 1 ratio, and the other two have 1 : 2 ratio of the 
number of positively charged atoms to the number of negatively 
charged atoms. Explicitly, the three typical types of surfaces of 1 : 1 
ratio are the type 2 surface obtained from (100) surface of a face- 
centered cube (such as the NaCl(100) surface), the type 3 surface 
obtained from (110) surface of simple cube (such as the CsCl(llO) 
surface) and the type 4 surface obtained from (110) surface of face- 
center cube (such as the NaCl(llO) surface), which are shown in 
Fig. 4(a) (b)(c), respectively. The other two typical types of surfaces 
of 1 : 2 ratio are the type 5 surface obtained from (110) surface of a 
face-centered cube (such as the CaF 2 (110) surface) and the type 6 
surface obtained from (110) surface of a body centered tetragonal 
(such as the BaO 2 (110) surface), as shown in Fig. 4(d)(e), respect- 
ively. For the two typical types of surfaces of 1 : 2 ratio, in order to 
keep the system charge neutral, we set the quantity of each positive 
charge twice of the quantity of each negative charge. Thus, for easy 
description, the notion q in the contexts below corresponds to the 
quantity of the positive charge for the two types of 1 : 2 ratio. The 
structure shown in Fig 1(a) is denoted as type 1. In Fig. 1(d), we show 
the critical lengths of the charge dipoles for various types of the solid 
surfaces. One can observe that the explicit value of the critical value 
depends on the geometry of the solid surface and the dipole struc- 
tures. Explicitly, the critical length of the charge dipoles on the solid 
surface is l c = 0.202 nm for q = 1.0 e and l c 0A = 0.302 nm for q = 
0.4 e for the types 1-3. For the type 4, the critical length of the charge 
dipoles on the solid surface is l c — 0.162 nm for q = 1.0 e and l c 0A = 
0.262 nm for q = 0.4 e, which are smaller than the case of type 1-3. 
As for the type 5 and 6, the critical length of the charge dipoles on the 
solid surface is / C L0 = 0.162 nm for q = 1.0 e and l c 0A = 0.302 nm for 
q = 0.4 e; Z C L0 = 0.202 nm for q = 1.0 e and l c 0A = 0.362 nm for q = 
0.4 e, respectively. The close values of the critical lengths in the solid 
surfaces of the first three types may be attributed to that there is only 
one characteristic length for each type. For the type 4 solid surface, 




Figure 4 | Other five typical surfaces of binary compound crystals, which 
correspond to the surface type number = 2, 3, 4, 5 and 6, respectively. The 
green and blue spheres represent positively and negatively charged atoms, 
respectively. 



there are two characteristic lengths, where the longer one is \l times 
larger than the shorter one. Here we only consider the shorter one as 
the critical length so that the critical length shows a smaller value. As 
for the types 5 and 6, the number of positively charged and negatively 
charged atoms is different and there are also two characteristic 
lengths so that the critical lengths show different behavior from those 
in the surfaces of the first three types. The lengths of the chemical 
bonds of some typical binary compound crystals or chemical groups 
are shown in Fig. 1(d), including the lengths of the C-O bond in 
-COOH group, C-N bond in -CONH 2 group, and the bonds of 
BN, CuO, A1 2 0 3 , CaF 2 and CaO. It is clear that all of them are less 
than l c ° A and most of them are less than l c . This demonstrates that 
the observations should be of general importance. 

In summary, we have shown that the wetting properties of a solid 
surface significantly depend on the surface topology parameter, 
which can be characterized by the length of the charge dipoles /. 
Counter to intuition, the length of the charge dipoles on the solid 
surface plays a key role on the wetting behavior and there is a critical 
length of the charge dipoles, below which the charge dipoles on the 
solid surface have a negligible role in the wetting property. The solid 
surface still exhibits hydrophobic behavior when the dipole length is 
less than the critical value, indicating that the water molecules on the 
solid surface seem not to "feel" attractive interactions from the 
charge dipoles on the solid surface, no matter how large the moment 
magnitude of the charge dipoles is. Those unexpected observations 
result from the collective interactions between the water molecules 
and charge dipoles on the solid surface, where the steric exclusion 
effect between water molecules prevents those hydrogen atoms of 
water molecules from staying very close to the negative charge and 
those oxygen atoms of water molecules from staying very close to the 
positive charge, reducing the interactions between the water mole- 
cules and the charge dipoles. Remarkably, the steric exclusion effect is 
also important for surfaces with charge dipole lengths / greater than 
this critical length. When charge dipole lengths are greater than this 
critical length, the surfaces become more hydrophilic with the dipole 
lengths increase. We note that the bond lengths of most of common 
materials are less than the critical lengths, suggesting that the obser- 
vation should be of general importance to the surface wetting prop- 
erties. We further note that the critical length of the surface still exists 
for different Lennard- Jones parameters of the atoms on the substrate 
within a reasonable range (see Table S2 in the Supplementary 
Information). Moreover, when the charge dipoles are provided from 
planting groups, such as -COOH, the critical values will be smaller 
due to the flexible behavior of the groups, which decreases the steric 
exclusion effect. These findings also have biological significance 
and are helpful for interfacial engineering, new functional material 
design including carbon and silicon nanostructures, and offer 
insights into many unexpected observations on wetting behavior. 

Methods 

A constant temperature and constant volume (NVT ensemble) MD simulation was 
performed using a time step of 1.0 fs with Gromacs 3.3. 1 58 . The periodic boundary 
conditions were applied in all directions. The selection of a vapor-liquid coexistence 
system was used to maintain the ambient conditions, and a Berendsen thermostat 59 
with a time constant of 1.0 ps for coupling was used to maintain the temperature of 
the water at 300 K. The solid atom with Lennard- Jones parameters e ss = 0.105 kcal/ 
mol and o ss = 3.343 A, and the SPC/E water model was used as the explicit solvent. A 
cutoff of 10 A was used for both the particle-mesh Ewald method 60 with a real space to 
model long-range electrostatic interactions, and for the van der Waals interactions. 
The contact angle of a water droplet was calculated from the density profile. For 
simplicity, when all of the water molecules completely spread on the surface, we 
assumed that the contact angle was 0. In the computation of the radial distribution 
function (RDF) of the water molecules next to a given surface and the interaction 
between those water molecules and the surface, we considered the solid surface with 
1664 atoms immersed in the center of a water box with 4 nm thickness of water layer, 
where periodic boundary conditions apply along the three dimensions 52 . 
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